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ABSTRACT 

Aims. We aim at obtaining very accurate positional constraints for seven gravitationally lensed quasars currently monitored by the 
COSMOGRAIL collaboration, and shape parameters for the light distribution of the lensing galaxy. We also want to find simple mass 
models that reproduce the observed configuration and predict time delays. Finally we want to test, for the quads, whether there are 
clues of astrometric perturbations coming from substructures in the lensing galaxy, preventing from finding a good fit of the simple 
models. 

Methods. We apply the iterative MCS deconvolution method to near-IR HST archives data of seven gravitationally lensed quasars. 
This deconvolution method allows us to separate the contribution from the point sources to the one from extended structures such as 
Einstein rings. That method leads to an accuracy of 1-2 mas on the relative positions of the sources and lens. The limiting factor of 
the method is the uncertainty on the instrumental geometric distortions. We then compute mass models of the lensing galaxy using 
state-of-the-art modeling techniques. 

Results. We obtain relative positions for the lensed images and lens shape parameters of seven lensed quasars: HE 0047-1756, 
RXJ1131-1231,SDSSJI138+0314, SDSS JI 155+6346, SDSSJ1226-0006, WEI J2026-4536 and HS 2209-1-1914. The lensed image 
positions are derived with 1-2 mas accuracy. Isothermal and de Vaucouleurs mass models are calculated for the whole sample. The 
effect of the lens environment on the lens mass models is taken into account with a shear term. Doubly imaged quasars are equally 
well fitted by each of these models. A large amount of shear is necessary to reproduce SDSS Jl 155+6346 and SDSS J1226-006. In 
the latter case, we identify a nearby galaxy as the dominant source of shear The quadruply imaged quasar SDSS Jl 138+03 14 is well 
reproduced by simple lens models, which is not the case for the two other quads, RX Jl 131-1231 and WEI J2026-4536. This might be 
the signature of astrometric perturbations due to massive substructures in the galaxy unaccounted for by the models. Other possible 
explanations are also presented. 

Key words. Gravitational lensing: astrometry - Quasars: general - Cosmology: observations - Techniques: image processing 

1 . Introduction lution frames. A good example of the consequences on Hq of ac- 

I I . . curate astrometry is t he quadruply lensed quasar PG 1 1 15+080: 

iReftdaH ilSM) was the first to state that gravitationally lensed ICourbin et alJ (fl997h found stronger constraints on the lensed 

quasars can be very useful for determining parameters of our j^^ges and lensing galaxy positions which reduced the degen- 

Universe: combined with a model of the mass distribution in the grades and the range of acceptable models. This allowed to mea- 

lensing galaxy, the time delay between different lensed images 5^.5 //„ with an accuracy two times better than in previous stud- 

H can lead to the determination of the Hubble constant, Hq. This je^. fKeeton & KochanelJ (fT997h also showed that reducing the 

. . = motivated many of the early lensed quasars studies and time de- error on the lens galaxy position of PG 1 1 15+080 from 50 mas 

lay measurements campaigns. Unfortunately it quickly became jq 10 mas could impr ove the constraints on the lens model and 

cleai- that the final estimate of Ho is very sensitive to systematic j^us on Hq. Moreover lLd^^TTZI dSl) highlighted that a poor 

errors on the lens modeling. One way to reduce these systematic knowledge in the position of the lens galaxy in two different sys- 

errors is to derive accurate relative astrometry of gravitationally jgms, i.e. B0218+357 and PKS 1830-21 1, prevented from accu- 

lensed images and lens galaxy light profiles based on high reso- ^.^^^^y measuring the Hubble constant. Systematic studies of the 

effect of the astrometric accuracies on Hq are difficicult as they 

* Based on observations made with the NASA/ESA HST Hubble depend on the lens system configuration. 
Space Telescope, obtained from the data archive at the Space Science 

Instimte, which is operated by AURA, the Association of Universities Strong lensing is also a promising tool to estimate the 

for Research in Astronomy, Inc., under NASA contract NAS-5-26555. amount (possibly as a function of redshift) of dark matter clumps 

** Research Eellow, Belgian National Eund for Scientific Research (hereafter, following other authors, we will call them "sub- 

(ENRS) structures") in distant galaxies an d compare it to predictions 

Alexander von Humboldt Eellow of numerical simulations (see e.g. IZackrisson & Riehml 120091 : 
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iKoopmans et al.ll2009ah . The first evidence that strongly lensed 
quasars are sensitive to substructures in galaxies comes from the 
so-called "anomalous flux ratios"; for many systems, the flux ra- 
tios between the lens ed images deviat e from those predicted by 



tios between tne lens ea unages aeviat e irom tnose preaicteu by 
simple lens models (lKochaneklll991l: [Mao & SchneideHll998t 



iDalal & Kochanekll2002t iKeeton et al.ll2003h . It has long been 
thought that substructures only act on the image flux ratios be- 
cause of the dependence of the latters on the second deriva- 
tive of the gravitational potential. However, recent works ex- 
plore two new routes to detect substructures in lensed quasars. 
On e method suggests to u s e time delay measurements, as shown 
by iKeeton & MoustakasI (l2009l) who have highlighted small 
changes in time delays because of substructures. Even if these 
delays are likely to be modified only by a few tenths of a per- 
cent, future large monitoring campaign s should allow the detec - 
tion of the signature of substructures (iMoustakas et al.|[2009h . 
Another method proposes to detect substructures in the lens- 
ing galaxy through their effects on the position of lensed im- 
ages. The amplitude and probability we should expect for this 
phenomenon is still debated. On one hand, observable astro- 
metric perturbations should be due to the most massive sub- 
structures. But becaus e of the scarci ty of high mass dark mat- 
ter clumps, .Metcalf & Madaul (1200 lb derived a low probability 
and, on average, low a strometric perturbations. On the other 
hand, IChen et aP (l2007 l) showed that lower mass substructures 
also play a role. Including a large range of sub-halo masses, 
they find that substructures c ould induce astrometric per turba- 
tions as large as 10 mas (see IZackrisson & Riehml 12009 1 for a 
more complete review). Observationaly, astrometric perturba- 
tions caused by substructures were detected in a few systems, 
the most remarkable ones being MG2016+ 112 ( Koor)man set al,| 
[2002; More et al. 2009) and B0128-I-437 (Biggs et al. 2001). In 
both cases, the anomalies have been unveiled thanks to high res- 
olution radio images. 

Although we cannot yet reach the spatial resolution of the 
Yery Long Baseline Array (V LBA) in the optica l rang e, it 
has been shown by some of us (IChantrv & Masai ni l2007h that 
a sophisticated deconvolution technique (ISMCSOl) applied to 
Hubble Space Telescope (HST) images could lead to relative po- 
sitions of lensed quasar images with milliarcsecond (mas) accu- 
racy, reducing the error bars by a factor > 2 compared to other 
techniques. In the present paper, we apply this technique to a 
sample of seven gravitational lenses without measured time de- 
lays. All these systems are photometrically monitored by the 
COSMOGRAIL0 collaboration and should get a time delay mea- 
surement in the near future. The goals of this paper are twofold. 
First we want to provide shape parameters for the lensing galaxy 
and accurate relative astrometry for these systems together with 
simple lens models. From the latter, prospective time delays 
are also calculated, comple menting time del ays predicted with 
non parametric modeling (iSaha et al.l l2006h . Second, we sys- 
tematically investigate, for quadruply imaged quasars, the ability 
of simple smooth models to reproduce the image configuration 
within a few milliarcsec. From this systematic and uniform ap- 
proach, we want to test whether the actual data show evidence 
for astrometric perturbations due to substructures. 

The lens sample studied in this paper is composed of 7 dif- 
ferent systems without time delay measurements and currently 
monitored by the COSMOGRAIL collaboration: 4 doubly im- 
aged quasars for which no detailed modeling and/or relative 



astrometry has ever been published and 3 quadruply imaged 
quasars, amongst which 2 have not yet been studied in detail 
either (no modeling, no time delay and/or no lens redshift). The 
ISMCS deconvolution of the gravitational lenses for which time 
delays have already been measured will be presented in another 
paper (Chantry et al., in preparation). 

The studied sample is detailed in Sect.[2]and the data in Sect. 
[3] while the image processing technique is explained in Sect. 
ID along with the results. The modeling strategy is explained in 
Sect. |5] A discussion of the models is presented in Sect. |6l We 
then conclude in Sect.|2] 



2. An overview of our sample 

Here are the seven gravitationally lensed quasars of our sample, 
the right ascension and declination being expressed in the J2000 
coordinates system: 

- HE 0047-1756 (a) 

This object (RA ^ 00''50'"27^82 and DEC = -17°40'08'.'79) 
was discovered bv lWisotzki et al.l (2000) in the framework of 
the Hamburg/ESO Survey (HES) for bright quasars, cover- 
ing the Southern sky. It was later identified by Wisotz ki et al.l 
(2004) as a doubly imaged quasar at a redshift of Zj = 1.68. 
The lens is an elliptical galaxy with a spectroscopic redshift 
of zi = 0.407 + O.OOI dEigenbrod et"al]|2006bt lOfek et al.1 
I2OO6I) . 

- RXJ1131-1231(b) 

This quadruply imaged quasar (RA - ll*'31™55i'.39 
and DEC = -12°31'54'.'99) was discovered in l2003l by 
ISluse et al. 'Thev found a redshift of zi = 0.295 ±0.002 for the 
lens while the source lies at Zj = 0.657 + 0.001. Preli minary 
time delays have been proposed by Morgan et al. (l2006h 
and revised est imates will be published in .Kozlovski et al.l 
din preparationh . The system was charac terized in d e tails i n 
terms of astrometry and photometry bv ISluse et alJ (l2006l) . 
Claeskens et al. (2006) modeled it and also reconstructed the 
source which appears to be a Typ e 1 Seyfert spiral galaxy . 
A similar work was perfomed by iBrewer & Lewis! d2008h 
using a Bayesian approach. Evi dence for substruct ures in 
the main l e ns wa s se arched by ISugai et aP (l2007h while 
ISluse et al.1 d2007h and iDai et all d2010h used microlensing 
to study the quasar source. 



' Iterative Strategy combined with the MCS deconvolution alogrithm 
^ COSmological MOnitoring of GRAvItational Lenses; 



|http : //www ■ CO smogr ai 1 ■ or g 1 



SDSSJl 138+0314 (c) 

This quadruply imaged object (RA = llh38'"03^70 and 
DEC = +03°14'57'.'99) was dis covered in i2008 in the Sloan 
Digital Sky Survey (SDSS) bv llnada et al.1 The redshifts of 
the quasar and the lens were measured by lEigenbrod et al.l 
(2006b) and are respectively equal to Zj = 2.438 and 
zi - 0.445 + 0.001. No detailed modeling and no time delay 
have ever been published for this system. 

SDSSJ1155+6346(d) 

This doubly imaged quasar (RA = ll'' 55"'I7^34 and DEC 
= -h63°46'22'.'00) was discovered by iPindor et all d2004 
in the SDSS data set. They measured the redshifts of the 
quasar and the lens: Zj - 2.888 and zi = 0.176. They also 
found that one of the two images of the quasar is very 
close to the lensing galaxy (at around 10% in effective 
radius off the center of the lens) and is the brightest. That 
configuration cannot be reproduced by a simple model of 
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mass distribution. 

SDSS J1226-0006 (e) 

This system (RA = 12''26'"08^02 and DEC = 
-00°06'02'.'19) is a doubly ima ged quasar disco vered 
in the framework of the SDSS by llnada et alj (|2008|) . The 
quasar is located at a redshift of Zj = 1.125. According to 
lEigenbrod et al. (2006b), the lens is likely to be an early-type 
galaxy, with a spectroscopic redshift of z\ = 0.516 + 0.001. 
This system has no measured time delay, no published 
relative astrometry and no detailed modeling study. 

WFIJ2026-4 5 36(f) 

[Morgan et al.l (|2004|) discovered this quadruply imaged 
quasar (RA = 20''26™10^43 and DEC = -45°36'27'.'10) 
during an optical survey using the WFI camera mounted on 
the MPG/ESO 2.2m telescope operated by the European 
Southern Observatory (ESO). The redshift of the source is 
Zs - 2.23. The one of the lens is unknown, although it is 
clearly visually detected on high resolution images. No time 
delay h as ever been measured but according to Morg anet al.l 
(l2004l) . the longest one might be of the order of at most 
a week or two. According to them, the lensed images are 
likely affected by microlensing. 

HS 2209+1914 (g) 

This system (RA = 22*^1 1'"30^30 and DEC = 
+ 19°29'12'.'00) is a doubly imaged quasar, with z, = 1.07, 
discovere d during the Hambu rg-Cfa Bright Quasar Survey 
(HS) by Hag en et akl (11999) . They clearly detected the 
lensing galaxy. Nothing else is available for this system: no 
time delay, no lens redshift and no modeling. 



3. Observational material 

The images we analyse were acquired with the camera 2 
of NICMOS, i.e. the Near-Infrared Camera and Multi-Object 
Spectrometer (hereafter NIC2) mounted on the HST. They were 
all obtained in the framework of the CASTLES project (Cfa- 
Arizona Space Telescope LEns SurvejQ, PI: C.S. Kochanek), 
and are available in the HST archives. The filter used is the 
F160W which is very close to the //-Band filter It was selected 
for several reasons: on one hand, the PSF is well-sampled and 
only slightly variable across the field; on the other hand, the 
lensed quasars flux ratios are less aff'ected by microlensing ef- 
fects, given the size of the quasar at these wavelengths, and dust 
extinction by the lensing galaxy. Details about the image acqui- 
sition are summarized in the first columns of Table [T] the name 
of the object, the date of observation, the number of frames 
and the total exposure time. All the frames were obtained af- 
ter the installation of the NICMOS Cooling System, or NCS, 
in 2002. Every image was acquired with dithering and in the 
MULTIACCUM mode, each one of them being a combination 
of about twenty subframes. As these objects were all observed 
between October and december 2003, the pixel size of the de- 
tector on the sky does not change from one target to the other, 
also because the plate scale of NICMOS has become very stable 
since the installation of the NCS. The values we use were mea- 
sured during part b of the third Servicing Mission Observatory 
Verification, SM0V3b, and are the foUowiri g: x=0.075948"and 
v^0.075355" (STScI NICMOS GroudlloOTi) . 



4. ISMCS on HST/NIC2 images 

4.1. Technique 

To extract accurate spatial and shape parameters from our data, 
we need a method capable of separating the contributions of the 
lensed point sources from the ones of the more diffuse compo- 
nents (galaxies, halos, arcs, ring s, ...). This is exactl y what the 
MCS deconvolution algorithm (iMagain et al.l 1 19981) provides. 
One of the advantages of this deconvolution method with respect 
to other techniques is that it does not violate the sampling the- 
orem. In practice that means that we do not try to fully decon- 
volve an image in order to obtain an infinite resolution. Instead, 
we choose a resolution for the final deconvolved image, in 
our case a Gaussian with 2 pixels Full-Width-at-Half-Maximum 
(FWHM), and we deconvolve our images with a partial Point 
Spread Function or PSF (which gives the total PSF when re- 
convolved with our 2 pixels FWHM Gaussian). To achieve this 
task, we need to know very well the shape of the PSF. As the 
NIC2 field is only 19'.'2 x 19'.'2, we do not have the possibil- 
ity to use field stars to determine the PSF. Moreover, since the 
lensed quasar images are contaminated by the lensing galaxy or 
partial Einstein rings underneath them, we cannot use these im- 
ages direc tly to improve our PSF. Instead, we use ISMCS (see 
IChantrv & Magain 2007, for further details), a special iterative 
strategy coupled with the MCS algorithm. The HST PSF being 
quite complex (it includes spike-like features and an intense first 
Airy ring), we start the deconvol ution process using a PSF cre- 
ated by the Tiny Tim software ( "Krist & Hookl 12004 !) as a first 
guess of the true PSF. We improve the Tiny Tim function in ad- 
justing it at the same time on all the point s ources of a frame us- 
ing a technique described in iMagain et al.l (.2007.) which allows 
to add a numerical component to the input PSF so that it is bet- 
ter adapted to the actual frame. We then obtain a set of modified 
PSFs which we use to simultaneously deconvolve all the frames 
at our disposal, with a sampling step two times smaller than the 
original one. In doing so, we obtain a first approximation of the 
diffuse background and after reconvolving it to the inital resolu- 
tion, we subtract it from the original frames and we obtain new 
ones, partially cleaned from the extended structures. On these 
modified frames which contain point sources less contaminated 
by smooth structures, we improve once again our PSFs. This it- 
erative process has to be repeated until the reduced chi square, 
X^ (Eq. [TJ, reaches a value close to unity in an area determined 
by the maximum extension of the residual structures after the 
first deconvolution, and until the residuals are sufficiently flat 
(no sharp structure). In practice, we stop when an additional it- 
eration no longer improves the;i'^ (typically, A^^ < 0.2). For an 
image with A^ pixels, the latter is defined as follows: 



^^ N lu 



M{x) - D(x) 
o-(x) 



(1) 



|http : //www ■ cf a . harvard . edu/castles/| 



where Ai(x) is the model reconvolved by the partial PSF, !D(x) 
is the observed signal and cr(x) the standard deviation associ- 
ated to that signal. In practice and as an improvement compared 
to the process applied to the Cloverleaf in the original paper 
of Chantrv & Magain (2007), we noticed that convergence is 
reached faster when performing a first simultaneous deconvolu- 
tion of all the frames with Tiny Tim PSFs instead of first trying 
to improve these PSFs on the unmodified images, i.e. containing 
the untouched background structures. With this first deconvolu- 
tion, we obtain a map of the diffuse structures. The latter has then 
to be cleaned from some artificial ring structures present only to 
compensate for inaccurate PSFs. It can then be subtracted from 
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Object 


Date of obs. (y-m-d) 


# of frames 


Total exp. time 


# of iterations 


;^f^ first it. 


x'-j last it. 


(a) HE 0047-1756 


2003-12-10 


4 


44' 


3 


60.49 


1.39 


(b)RXJ1131-1231 


2003-11-17 


8 


89' 


4 


60.18 


2.26 


(c)SDSSJl 138+0314 


2003-11-06 


4 


44' 


3 


19.61 


1.76 


(d)SDSSJl 155+6346 


2003-12-12 


5 


84' 


4 


30.15 


3.33 


(e)SDSSJ1226-0006 


2003-11-21 


4 


44' 


3 


20.03 


1.08 


(f) WFI J2026-4536 


2003-10-21 


4 


46' 


4 


36.62 


4.35 


(g)HS 2209+1914 


2003-10-14 


4 


44' 


4 


40.39 


4.2 



Table 1. General information about the acquisition of HST/NIC2 images and the apphcation of ISMCS. 



the original frames which can then be used to improve a first 
time the PSFs. That first deconvolution is accounted for as one 
iteration. 



4.2. Results 

The original frame (combination of all observations), the decon- 
volved imageQ obtained from the last iteration of ISMCS and 
the mean residual mapfl after the first and last iterations are dis- 
played for each system in Fig. [T] Both residual maps are ex- 
pressed in units of cr and their color scale ranges from -5 in 
black to +5 in white. The black rectangle delimits the zone taken 
into account to estimate the reduced ;t^^, the orientation of these 
rectangles being the one of the original frames. We emphasize 
that the PSF used for the first iteration is the one created by the 
Tiny Tim software. When we examine the residual maps, the im- 
provement brought by ISMCS is undeniable. Moreover, in most 
cases (5 amongst 7), the remnant structures underneath the point 
sources (on the residual map from the last iteration), are in dis- 
agreement with each other, which is the sign of a variable PSF 
throughout the detector, even on small spatial scales. The num- 
ber of iterations necessary to reach convergence is shown in the 
last columns of Table [T] along with the values of the x^ after the 
first and the last iterations. 

The astrometry, corrected from the X/Y scale difference and 
the distortions of NIC2, and the photometry (Vega system) are 
shown in Table |2] The +lcr error bars were calculated in decon- 
volving each frame individually at the last iteration and in de- 
termining the dispersion around the mean. They are very small 
because they are inherent to the deconvolution technique: no ex- 
ternal systematic error is included in these error bars. To estimate 
the total error, we compare the spatial extension of each object 
on the detector to the one of the Clo verleaf (H 14 1 3+ 1 17). Th e 
latter was used as a test of ISMCS in Chantry & MagainI (l2007h : 
in comparing the astrometry of the point sources obtained in two 
different filters and with two different orientations on the sky, 
they could estimate the total error to 1 mas, accounting e.g. for a 
possible remnant distortion in the images. The estimated total er- 
rors based on the Cloverleaf are displayed in the fifth column of 
Table|2l Of course, as they are based on the maximum extension 
of the object no matter the direction, they should be considered 
as upper limits. That is why they are called "MTE" which stands 
for "Maximum Total Error". 

Since the total error derived in lChantrv & MagainI ilOOl) for 
H1413+1 17 was based on a comparison of the relative positions 
of the lensed images obtained at different NIR wavelengths and 
image orientations with the same instrument, we have attempted 



to get independent estimates based on the comparison of the 
relative astrometry derived with HST and with high resolution 
radio data. In a future paper treating the lenses with akeady 
measured time delays (Chantry et al., in preparation), we will 
present HST astrometry for the radio quad JVAS B 1422+231 
iPatnaik et al. 1992). To estimate the error affecting our results, 
we choose one lensed image as astrometric reference and we 
calculate the distance between it and every other lensed image. 
We then measure the difference between the distances obtained 
with our posit ions an d the ones calculated with the radio astrom- 
etry of Patnai ket alJ (11999) . The scatter of these differences of 
distance around the mean is about 2.6 mas. Assuming the un- 
certainty is the same in any direction, we derive an error on the 
relative astrometry of 1.8 mas in RA and in DEC. This value is 
larger then the MTE of 1.05 mas derived for B 1422+231 from 
our standard method. This is expected as the radio emission in 
B 1422+231 is slightly extended and is likely not to originate 
from the accretion disk (observed in the optical range) but rather 
from a nearby region at the basis of the radio jet. Such an effect, 
known as core-shift, is observed betwe en two different radio- 
bands (lPorcasll2009HKovalev et alj|2008l) and may induce astro- 
metric perturbations as large as a few mas on the r elative as- 
trometry of lensed quasar images ( Mittal et alj|2006 ). Thus, it 
appears that the use of radio data as an independant calibrator of 
the systematic errors at mas accuracy is difficult. It requires the 
comparison of the relative astrometry, for several objects, be- 
tween radio and optical wavelengths, a task that is beyond the 
scope of the present paper 

Together with the point-source deconvolution, we use an an- 
alytical model to characterize the lensing galaxy light distribu- 
tion. To ensure in this case that the maximum amount of light 
of the galaxy is included in the profile, the deconvolution is per- 
formed with no numerical component. Since most of the lens- 
ing galaxies are ellipticals, we use a de Vaucouleurs light pro- 
file (de Vaucouleurs 1948). This procedure allows us to extract 
the galaxy shape parameters summarized in Table [3j the posi- 
tion angle or "PA" (orientation in degrees East of North) of the 
galaxy, its ellipticity (e = 1 - b/a), the effective semi-major and 
semi-minor axis (resp. Ocff and befj). The effective radius R^ff 
is further calculated as b eing the geometri cal mean between the 
two effective semi-axis JKochanekl l2002h . These three quanti- 
ties are expressed in arcseconds. The ±lcr error bars were also 
calculated in deconvolving each frame individually and in de- 
termining the dispersion around the mean. Let us note that the 
luminosity of the galaxies displayed in Table|2]is measured in an 
aperture equal to Reff- 



'' The labels of the lensed images are the same as in previous studies 4.3. DiSCUSSIon 
if any. 

^ The residual map is the image of the difference between the model A few remarks can be made about the results from the deconvo- 

and the original frame in units of sigma. lution: 



V. Chantry et al.: COSMOGRAIL: Accurate astrometry and models for 7 lensed quasars 



Object 


Label 


ARA(") 


ADEC (") 


MTE (mas) 


Magnitude 


Flux ratio 


(a) HE 0047-1756 


A 


0. 


0. 


1.17 


15.19 ±0.01 


1. 




B 


0.2328 ± 0.0008 


-1.4094 ±0.0002 


1.17 


16.69 ±0.01 


0.253 ± 0.002 




G 


0.2390 ± 0.0022 


-0.8098 ± 0.0056 


1.17 


17.17 ±0.02 


/ 


(b)RXJl 131-1231 


A 


0. 


0. 


2.64 


15.36 ±0.01 


1. 




B 


0.0347 ± 0.0005 


1.1870 ±0.0005 


2.64 


15.58 ±0.01 


0.816 ±0.003 




C 


-0.5920 + 0.0007 


-1.1146 ±0.0004 


2.64 


16.42 ± 0.01 


0.374 ± 0.003 




D 


-3.1154 ±0.0012 


0.8801 ± 0.0013 


2.64 


17.76 ±0.01 


0.110 ±0.001 




G 


-2.0269 ± 0.0016 


0.6095 ± 0.0015 


2.64 


15.55 ±0.03 


/ 


(c)SDSSJ1138+0314 


A 


0. 


0. 


1.17 


17.89 ±0.01 


1. 




B 


-0.1003 ±0.0006 


0.9777 ± 0.0007 


1.17 


19.07 ± 0.01 


0.336 ± 0.004 




C 


-1.1791 ±0.0003 


0.8119 ±0.0007 


1.17 


18.89 ±0.01 


0.400 ± 0.002 




D 


-0.6959 ± 0.0003 


-0.0551 ±0.0003 


1.17 


19.02 ± 0.01 


0.354 ± 0.003 




G 


-0.4633 ± 0.0071 


0.5340 ± 0.0036 


1.17 


17.77 ±0.01 


/ 


(d)SDSSJl 155+6346 


A 


0. 


0. 


1.59 


16.83 ± 0.02 


1. 




B 


1.8983 ±0.0005 


0.4052 ± 0.0005 


1.59 


17.87 ±0.01 


0.710 ±0.017 




G 


1.6982 ±0.0024 


0.3438 ± 0.0009 


1.59 


15.71 ±0.01 


/ 


(e) SDSS J1226-0006 


A 


0. 


0. 


1.03 


17.05 ± 0.01 


1. 




B 


1.2563 ±0.0002 


-0.0550 ± 0.0007 


1.03 


17.80 ±0.01 


0.499 ± 0.006 




G 


0.4386 ± 0.0029 


0.0209 ± 0.0034 


1.03 


17.71 ±0.03 


/ 


(f)WFIJ2026-4536 


B 


0. 


0. 


1.17 


17.08 ±0.01 


1. 




Al 


0.1613 ±0.0007 


-1.4290 ±0.0005 


1.17 


15.58 ±0.01 


3.988 ±0.018 




A2 


0.4140 ± 0.0007 


-1.2146 ±0.0006 


1.17 


16.03 ±0.01 


2.634 ± 0.017 




C 


-0.5721 ± 0.0006 


-1.0437 ±0.0003 


1.17 


17.26 ± 0.01 


0.851 ±0.07 




G 


-0.0479 ± 0.0015 


-0.7916 ± 0.0015 


1.17 


18.94 ±0.04 


/ 


(g)HS 2209+1914 


A 


0. 


0. 


0.85 


14.37 ± 0.02 


1. 




B 


0.3307 ± 0.0004 


-0.9863 ± 0.0010 


0.85 


14.63 ±0.01 


0.790 ± 0.027 




G 


0.2155 ± 0.0037 


-0.3947 ± 0.0054 


0.85 


21.58 ±0.2 


/ 



Table 2. Relative position, maximum total error ("MTE"), magnitude and flux ratio of the lensed images and lensing galaxy (see 
Fig.[I]for the labels). 



Object 


PAD 


e 


a.ff (") 


b.ff (") 


Reff (") 


(a) HE 0047-1756 


113.8 + 5.5 


0.22 ± 0.02 


1.02 ± 0.03 


0.81 ±0.02 


0.91 ±0.02 


(b)RXJ1131-1231 


108.6 ± 2.4 


0.25 ± 0.04 


1.25 ± 0.06 


0.97 ± 0.01 


1.11 ±0.03 


(c) SDSS Jl 138+0314 


122.7 ± 6.5 


0.16 ±0.02 


0.93 ± 0.04 


0.79 ± 0.03 


0.86 ± 0.03 


(d) SDSS Jl 155+6346 


0.7 ± 3.4 


0.15 ±0.02 


1.23 ± 0.02 


1.06 ±0.01 


1.14 ±0.01 


(e) SDSS J1226-0006 


45.2 ±6.1 


0.07 ± 0.02 


0.72 ± 0.04 


0.67 ± 0.02 


0.69 ± 0.03 


(f) WEI J2026-4536 


60.8 ± 5.4 


0.24 ± 0.03 


0.72 ± 0.03 


0.57 ± 0.02 


0.64 ± 0.02 


(g)HS 2209+1914 


63.1 ±3.25 


0.05 ± 0.02 


0.55 ±0.01 


0.52 ± 0.01 


0.53 ±0.01 



Table 3. Measured shape parameters for the lensing galaxy. 



HE 0047-1756 (a) 

A faint Einstein ring, stretched image of the quasar host 
galaxy, is revealed by the deconvolution. 

RXJ1131-1231( b) 

ISluse et aP (l2006l) reports astrometric measurements on the 
same frames with the MCS deconvolution algorithm but 
with no iterative strategy. Their results agree within the 
error bars with those presented here. An offset of up to 3 
mas between both results is observed. This difference is 
probably due to the large brightness of the Einstein ring. 
Indeed, the dififerent amount of recovered background under 
the PSF can lead to a small shift in position. Photometry 
is also afffected by the presence of the ring. W e derive an 
absolu te photometry about 0.4 mag brighter than lSluse et al.l 
(l2006h but we obtain compatible flux ratios. The remnant 
systematic structures in the final residual map also result 
from the presence of this very bright ring which affects the 
PSFs and degrades their quality. Indeed, because a part of 
the background is identical under the three brightest lensed 
images (and thus with more weight in the determination of 



the PSF), it is impossible to completely disentangle the flux 
contribution of the ring from the one really included in the 
point sources. This is a limitation of the ISMCS method: for 
it to work properly, the background has to be different under 
each point source. 

SDSS Jl 138+0314 (c) 

A faint Einstein ring is revealed by the deconvolution 
process. 

SDSS J1155+6346 (d) 

Our astrometry is not in agreement with lPindoret al.l (l2004t) 
especially concerning the lens. The dilTerence for source 
B amounts to U.'l in RA and 0'.'04 in DEC while for the 
lens the offset is much larger and amounts to 1'.'55 in RA 
and 0'.'28 in DEC. However, our astrometry is in agreement 
with the one listed in the CASTLES database. The remnant 
systematic structures in the final residual map are due to 
the presence of the very bright and extended lensing galaxy 
(which is not clear on the presented frames, the cuts being 
chosen so that the two lensed images appear clearly). As in 
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the case of RX Jl 131-1231, it is not possible to completely 
disentangle the background flux contribution from the one 
of the point sources, a part of the background being identical 
under both lensed images. This degrades the quality of the 
PSFs and of the deconvolution. 

- WFI J2026-4536 (f) 

The astrometry we obtain, except for the right ascension of 
the lens (ARA =0'/ 03), is in agreement, within the error bars, 
with the results of Morgan et alJ (12004 ) who used the same 
frames but a different image processing technique. 

- HS 2209+1914 (g) 

A bulge is clearly observable but an additional extended 
structure is also visible. It could either be some spiral arms, 
in which case the lens would be a late-type galaxy, or even 
a distorted Einstein ring. A spectrum of the lens and higher 
resolution imaging would help to disentangle between these 
hypotheses. As we do not know what this structure is related 
to, we fit a de Vaucouleurs model on the bulge only, in using 
a special feature of the MCS algorithm: a mask encircling the 
lens galaxy to avoid the model to fit this extended structure. 
Moreover, the residual map contains many intense structures. 
However, as these structures do not have the same shape un- 
der both point sources, we cannot recover them with another 
iteration of ISMCS. Saturation is unlikely to be responsible 
for that phenomenon, as it is corrected by the NICMOS re- 
duction pipeline. It could thus be due to differential extinc- 
tion by the lensing galaxy, resulting in a different color for 
both lensed images and thus a different shape of the PSF. 

Let us note that all our results are in agreement with what can 
be found in the CASTLES database, within their error bars (ours 
being smaller). 



5. Parametric modeling 

Our goals are twofold. First, we aim at providing simple models 
and prospective time delays for the lensed quasars monitored by 
COSMOGRAIL. Second, we aim to test whether simple smooth 
lens models are able to reproduce the mas relative astrometry of 
quadruply imaged quasars in 3 systems without measured time 
delays. Our strategy c onsists in usin g the LENSMODEL soft- 
ware package vl.99o (lKeetonll200lh to model the mass distri- 
bution of our seven systems. For a chosen model of the mass 
distribution, the code minimizes a x^ defined as the square of 
the difference between observable quantities and their "model 
counterparts", weighted by the observational errors on these 
quantities. Two different lens models are considered. First, an 
isothermal profile, which i s the standard mass distr ibution to 
model gravitational lenses (iKassiola & Kovnenll993 '). and sec- 
ond, a de Vaucouleurs profile, for which we assume that the 
light perfectly traces the mass in the inner regions of lensing 
galaxies. These two models should provide a good approxima- 
tion of the extreme slopes of the mass distribution at t he location 
of the lensed images and of the expected time delays (iKochanekl 
[2002 ; Kocha nek & Schechter 2004) . In addition, the sfiidy of 
the galaxy-galax y lensing sample from the Sloan Lens ACSQ 
Survey (SLACS, Bolton et aT]l2006h has revealed that the mas- 
sive elliptical lensing galaxies are nearly kinem atically undistin- 
guishable from isothermal ellipsoids (see e.g. Koorima ns et al.l 
i2009b) . This supports the use of an isothermal gravitational po- 



ACS stands foi Advanced Camera for Surveys. 



tential as a fiducial model to test the ability of smooth lens mod- 
els to reproduce quadruply imaged quasars with mas accuracy. 
Since lensing galaxies are never isolated, we model the effect of 
the environment with an external shear term characterized by an 
amplitude y and a position angle 0y (pointing towards the mass 
at the origin of the shear). All the models are computed for a flat 
universe with the following cosmological parameters: Hq - 70 
km/s/Mpc, Q,„ = 0.3 and Q.^ = 0.7. 

To model our systems, we use every constraint at our dis- 
posal: the relative astrometry of the lensed images, with the 
MTE, i.e. the uncertainties displayed in the fifth column of Table 
|2] the position of the main lens, with the error inherent to the 
deconvolution when it is larger than the MTE, and, in the case 
of doubly imaged quasars, the flux ratio between the two point 
sources. In principle, fluxes can be contaminated by different 
effects such as microlensing by stars in the lens galaxy, dust 
extinction and also by the time delay itself. As the flux ratios 
are m easured in the near-i nfrared, all these effects should be 
small (lYonehara et alJl2008l) . We thus assume a Icr error of 10% 
on the flux ratios. In summary, we have 10 constraints for the 
quads, while we have 8 for the doubles. For the de Vaucouleurs 
model, we assume that the total mass profile follows the light 
profile. We thus add three constraints to the model: the PA of the 
galaxy, its ellipticity e and effective radius R^ff (see Table [3]l. 
Due to the limited number of observational constraints, isother- 
mal mass profiles are assumed to be spherically symmetric (SIS, 
i.e. Singular Isothermal Sphere) when modeling doubles. This is 
not a strong asumption as the quadrupole term of the potential 
modifies only slightly the time del ays of doubly imaged quasars 
('Kochanek'l2002VWucknitz!'2002). For quads, we allow the el- 
lipticity of the isothermal mass distribution (SIE, i.e. Singular 
Isothermal Ellipsoid) to deviate from the ellipticity of the light 
profile. This enables to account for dark matter halos that would 
be rounder then the light distribution (Ferreras et al. 2008). The 
position angle of the total mass distribution can be constrained 
as the one of the light profil e as these two distr i butions might 
only be slightly misahgned (iKeeton et al.l 119971: iFerreras et al.l 
12008"). Finally, we also assume that the center of the total mass 
distribution and the one of the light profile are identical within 
the error bars. This is supported by the work of Yoo et al. (20061) 
who found, for 4 lensed quasars with an Einstein ring, that the 
offset between the light and the total mass distribution is limited 
to a few mas. Calculating the number of degree(s) of freedom 
(d.o.f.), which is the difference between the number of model pa- 
rameters and observable quantities, we find d.o.f. when model- 
ing doubly imaged quasars and 2 (resp. 3) d.o.f. when modeling 
quads with SIE (resp. de Vaucouleurs) + external shear. 

The search for the best model and estimate of uncertainties 
is performed in two steps. First, we generate an initial sam- 
ple of 2000 different models with parameters distributed over 
the whole parameter space and optimize them. This method 
is efficient to find the best models and identify local minima. 
Then, in order to estimate the model uncertainties, we sam- 
ple the posterior probability distribution of the parameter space 
using an adaptive Metropolis Hastings Monte-Carlo Markov 
Chain (MCMC) algorithm. This technique is implemented in 
LENSMODEL and described in Fadely et al. (2010). In practice, 
an ensemble of 15 different chains are run, each chain consisting 
of a sequence of trial steps drawn from a multivariate gaussian 
distribution of width estimated thanks to the first step of the pro- 
cess. The sampling of the parameter space is optimised by using 
the covariance matrix. In 5% of the steps, the covariance matrix 
is diagonal, allowing to use a large step along one of the axis 
and better escape local minima in the x^ surface. We use the 
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same criterion as lFadelvet alJ (1201 Ol) to assess that any MCMC 
run has converged. Finally, for each point of the MCMC, we 
calculate the relative likelihood of a parameter p based on the 
X^ statistics (i.e. L(D\p) = exp(-x^/2)), and calculate a 68% 
confidence interval for each parameter. 

The parameters of the best fit models are displayed in Table 
ID The columns display the following items: the name of the 
object, the type of mass distribution used ("DV" stands for 
de Vaucouleurs profile), the mass scale parameter (the angu- 
lar Einstein radius /?£,„ in arcseconds), the mass distribution 
ellipticity e and its orientation 9e in degrees positive East of 
North, the effective radius Reff in arcseconds in the case of a 
de Vaucouleurs model, the intensity of the shear y and its orien- 
tation 0y in degrees (East of North), the number of degree(s) of 
freedom (d.o.f.), the x^ of the fit and the predicted time delays 
in days when the lens redshift is known. For the quads and in the 
same column as i\\t x^^ we also give i\\tx\„ which is the contri- 
bution of the lensed images position to the x^i and x] which is 
the contribution of the lens galaxy position to the;^'^. Let us note 
that Mab > means that the flux of A varies before the one of B. 
The median value of each parameter along with 68% confidence 
level is shown in Table |5] 

6. Discussion 

6.1. Doubles 

For the doubly imaged quasars, both SlS-nshear and DV-nshear 
models can reproduce the image configuration as well as the 
flux ratio, even with our constraints on the shape of the galaxy 
in the case of a de Vaucouleurs profile. Two systems require 
a large amount of shear (i.e. y > 0.1 for both mass mod- 
els) to reproduce the lens configuration: SDSS J1226-0006 and 
SDSS J1155+6346. For SDSS J1226-006, the HST/N1C2 im- 
ages actually reveal a galaxy G2 at RA = r.'7153 and DEC - 
3'.'1710 from image A (3'.'4 from the main deflector), about 15° 
off the direction of 9y. This galaxy, which type is unknown, is 
likely not the only source of shear Indeed, the luminosity ratio 
between G2 and the lens is Lie„s/Lcn = 4.8. Assuming we ca n 
use the Faber-Jackson relation (L oc cr^. iFaber & Jacksonlll976h . 
this ratio leads to cri^„Jcr(j^ = 1.5, crig„s and 0"^, being respec- 
tively the velocity dispersion of the lens and of G2. The isother- 
mal model allows us to translate REin of the lens to crie„s. We 
find (Tie ns - 2\2 km/s and thus (Tg, = 141 km/s. Using formula 
A.20 of iMomcheva et all (l2006l) and supposing G2 is at the same 
redshift as the lens, this induces a shear of y = 0.039, more than 
2 times smaller than the one predicted by the SIS model. Other 
galaxies in the field are probably responsible for the remaining 
shear. 

A more dramatic case is SDSS Jl 155H-6346, for which mod- 
els predict a shear as large as 0.4 to reproduce the observed con- 
figuration. This is one of the largest shears needed to reproduce 
a lensed quasar system. On some larger field images of this ob- 
ject (obtained with ACS onboard HST, PI: C.S. Kochanek), we 
do not see any bright galaxy in its vicinity. We thus suspect that 
a massive galaxy cluster lies outside the ACS field, though noth- 
ing is clearly visible on the SDSS datcQ Deeper images would 
be necessary to infirm or confirm the existence of this cluster. 

In the case of HE 0047-175, a diffuse component lies at 
RA = -a.'0434 and DEC = -2'.'3393 from image A (1'.'56 from 
the lens), in the direction of the shear (see Fig. [1]) unregarding 
the employed model. Although very faint (about 2 mag fainter 
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than the lens), this galaxy is likely the major contribution to the 
shear in this system. Indeed, a SIS with cr = 88 km/s would pro- 
duce the observed amount of shear, if located at the position of 
this faint companion (assuming Zco^p - ziens = 0.407). 

6.2. Quads and astrometric anomalies 

The quadruply imaged quasars allow to test the ability of simple 
smooth models to reproduce a relative astrometry with mas accu- 
racy. Only the relative astrometry of SDSS J1138-H0314is easily 
reproduced with our models. Conversely, for WFI J2026-4536 
and RX J1131-1231, we find that a very large ;^^ is associated 
to our models. In the first case, the main contribution in the x^ 
comes from the difference between the RA^ of the model and the 
RA. of the light distribution. In the second case, the large x^ is 
mainly due to the impossibility of the model to recover the posi- 
tions of the lensed images and of the lensing galaxy. In any case, 
an underestimate of the error bars on the quasar lensed images is 
unlikely. 

For RX J1131-1231,a reduced ;(f^ ~ 1 can only be obtained 
if we increase the error bars on the positions of the lensed images 
by a factor 10. Alternatively, we also get a good fit if we allow 
more freedom to the position of the lensing galaxy (i.e. error = 
0.02). Following this procedure we find that the off'set between 
the light and mass distribution centroid amounts to 88 pc. This 
value is marginaly consi stent with the up per limit of 70 pc de- 
rived for B1938+666 bv lYooetall(l2006h . It is however incon- 
sistent with the maximum off'set value of ~ 20 pc found for the 3 
other systems they analysed, suggesting that the offset between 
light and mass distributions is not the caus e of th e astromet- 
ric perturbation we observe. IClaeskens et al.l (l2006l) also found 
that models (simple or more complex) fitting simultaneously the 
Einstein ring of RX J1 131-1231 and the quasar lensed images 
lead to a poorjt'-. Bre wer & Lew is (2008) were also unable to 
reproduce the lensed quasar relative positions to mas accuracy 
using the lens model based on the Einstein ring. 

For WFI J2026-4536, an acceptable x^ cannot be obtained 
in enlarging the error bars on the positions of the lensed im- 
ages or of the lens galaxy but only in relaxing the constraint on 
the position angle of this latter. This could be due to the pres- 
ence of a galaxy loc ated at RA = -7'.'398 and DEC = -1'.'940 
(iMorgan et alJl2004l) from image B, that we ignored in the mod- 
eling. 

For the three quads, the flux ratios predicted by our best mod- 
els differ from those measured on the HST/NIC2 frames. This 
discrepancy may have several origins such as invalid assump- 
tion about the lens model, i.e. need for multipole components 
(lEvans & Wit t^2003). micro lensing and/o r massive substructures 
in the lensing galaxies (see Keeton et alJ (12003. 2005) for an ex- 
haustive discussion). 

The amount of shear needed in our models to reproduce 
the configuration of the three quads is quite high (y > 0.1). 
Our model estimates can be compared to the minimum amount 
of shear required to reprodu ce the image positi on, following 
the methodol ogy described in Witt & Maol (Il997h . Using equa- 
tion (20a) of Witt & Mao ( Il997h . we find a minimum shear 
rniin of resp. 0.035, 0.004 and 0.062 for RX Jl 131-1231, 
SDSS Jl 138-1-0314 and WFI J2026-4536. The values found for 
the last two systems are much smaller than the one predicted 
by our models. This is not surprising as we find Oy and 0^ to 
be orthogonal, which implies that -y is in fact poorly known and 
strongly degenerate with the internal shear. For RX Jl 131-1231, 
Train is only 3 times smaller than the value derived for our SIE-t-y 
model. The difference between ymin and yobs is due to 6^ being 
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Object 


Model 


Rew 


e,e. 


^eff 


r,0r 


Flux ratios 


d.o.f. 


x' 


Time delays 


(a) HE 0047-1756 


SIS + 7 


0.751 


/ 


/ 


0.048, 7.36 


fB/fA = 0.253 





0.0 


AtAB = 11.9 




DV + 7 


0.756 


0.22, 113.78 


0.91 


0.120, 15.98 


fB/fA = 0.253 





0.0 


AtAB = 16.5 


(b)RXJ1131-1231 


SIE + 7 


1.834 


0.20, 117.52 


/ 


0.098, 96.37 


fB/fA = 0.615 
fc/fA = 0.553 
fo/fA = 0.053 


2 


200.5 
4„ = 66.0 
X] = 120.7 


AtAB = -l.0 

Mac = -1.3 
AtAD = 116.2 




DV + 7 


1.791 


0.32, 114.67 


1.10 


0.213, 101.73 


fB/fA = 0.679 
fB/fA = 0.584 
fB/fA = 0.042 


3 


184.6 
Xl = 56.9 

;r? = 118-4 


AtAB = -1.7 
AtAc = -2.4 
AtAD = 198.2 


(c)SDSSJl 138+0314 


SIE + 7 


0.6640 


0.05, 118.73 


/ 


0.107,32.12 


fB/fA = 0.505 
fc/fA = 0.714 
fo/fA = 0.945 


2 


2.5 

^L = o.i 

X^ = 1-9 


AtAB = 3.4 
AtAc = -1.7 
AtAD = 0.9 




DV + 7 


0.6629 


0.15, 121.40 


0.86 


0.145, 32.23 


fB/fA = 0.505 
fc/fA = 0.712 
fo/fA = 0.925 


3 


4.7 

^L = o.i 

X^ = 4.5 


AtAB = 3.8 
AtAc = -1.9 

AtAD = 1.0 


(d)SDSSJl 155+6346 


SIS + 7 


0.59 


/ 


/ 


0.392, 169.66 


fB/fA =0.710 





0.0 


AtAB = 20.6 




DV + 7 


0.58 


0.15,0.71 


1.14 


0.453, 168.98 


fB/fA =0.710 





0.0 


AtAB = 25.0 


(e) SDSS J1226-0006 


SIS + 7 


0.568 


/ 


/ 


0.100,8.01 


fB/fA = 0.499 





0.0 


AtAB = -25.5 




DV + 7 


0.557 


0.07,45.18 


0.69 


0.145,4.55 


fB/fA = 0.499 





0.0 


AtAB = -34.3 


(f) WFI J2026-4536 


SIE + 7 


0.6520 


0.22, 167.51 


/ 


0.151, 83.46 


,A,//b = 4.127 
/^J/s = 3.439 
/c//fi= 1.174 


2 


266.7 
Xl = 16.3 
X] = 66.2 


Zi unknown 




DV + 7 


0.6517 


0.25, 167.21 


0.64 


0.171,85.18 


/ai/./b = 3.998 
fAjfB = 3.307 
fc/fB = 1.096 


3 


263.1 
Xl = 14.4 
X^ = 62.9 


Zi unknown 


(g)HS 2209+1914 


SIS + 7 


0.515 


/ 


/ 


0.031, 94.27 


fB/fA = 0.790 





0.0 


Zi unknown 




DV + 7 


0.516 


0.05,63.10 


0.53 


0.041, 99.60 


fB/fA = 0.790 





0.0 


Zi unknown 



Table 4. Results of the parametric modeling. 



Object 



(a) HE 0047-1756 



(b)RXJ1131-1231 



(c)SDSSJ1138+0314 



(d) SDSS Jl 155+6346 



(e) SDSS J 1226-0006 



(f) WFI J2026-4536 



(g)HS 2209+1914 



Model 



SIS + 7 
DV + 7 



SIE + 7 
DV + 7 



SIE + 7 
DV + 7 



SIS + 7 
DV + 7 



SIS + 7 
DV + 7 



SIE + 7 
DV + 7 
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DV + 7 
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Table 5. Median value of the model parameters and 68% confidence interval. 



nearly aligned with 9y (about 15° offset), a situation which also 
leads to a sig nificant unde r estim ate of y using ymm- By using 
eq. (22) from IWitt & Mad (11997 ). we can also derive a range 
of allowed values for 9e based on the image configurations. For 
WFI J2026-4536, we find that the observed 9,, is far from the 
allowed values, confirming the results of our models which indi- 
cate a likely off'set betwen the mass and the light matter distribu- 
tion. For the two other systems, the observed 9^ falls at the limit 
of the allowed range, suggestive of a significant ellipticity of the 
lens, as observed. 

The previous results suggest that it is quite common for 
simple lens models to fail in reproducing mas astrometry of 
quadruply imaged quasars. To further investigate the question, 
we have searched the literature for lensed quasars having images 
with mas astrometric error bars (i.e. up to 0'.'002 on the lensed 
image positions) and published simple models. We found eleven 



( Keeton et al 



( Morga n etai] 120051: iKochanek etal. . 
\2006l lEigenbro detalT 



systems gathering these condition^ BO 128+ 437 dBiggs et al] 
200 |), MGq414+0534 dRos et al.j l2000l) . HE0435-1223 

l200^Sp SS0924+0219 
I2006ah H1413+117 
(Chantrv&Magaih 2007; MacLeod et al." ■2009"), B1422+231 
(Bradacetal. 2002). B 1 608+656 (Koopmans et al. 20031) 
B193 3+50^ (ICohnetal.1 [2001), MG20 16+112 dChen et al.1 
120071: iMore et al.' '200 91). WFI2033- 4723 dVuissoz et alJl2008l) . 
B2045+265 (McKean etal] l2007h . Two of these systems 
(H1413+117, B1933+503) are easily reproduced by simple 



** Although B0 712+4 72 has accurate astrometry, the model published 
by Jackso n et aU (11998) is provided without information about the x^ 
thus not allowing to estimate the quality of the fit. 

' B 1933+503 is actually a ten-images lens. Only the position of four 
of them is known with a precision of 2 mas or less. 
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models because of the large uncertainty affecting the position 
of the lensing galaxy (o-gai > (X'Ol). Out of the six systems for 
which VLA, VLBA or VLBI data are availabl^B (BO 128 +437, 
MG0414+0534, B1422+231, B1608-656, MG2016+112, 
B2045+265), only B1422+231 shows convincing evidence 
that smooth models allow to reproduce the relative astrome- 
try, although substructures are needed to reproduce the flux 
ratios (B radac et al.i i2002) . For MG0414+0534, B1608+656, 
MG2016+112 and B2045+265, complex models including a 
bright susbtructure (MG0414+0534) or a companion galaxy 
(B1608+656, MG2016+112, B2045+265) are needed to get 
to acceptable fits. The case of B0128+437 is a bit peculiar 
as the lensing galaxy was not detected at the time of the 
modeling pape£3 and so, only the lensed images were used 
for the models, the latter being constrained as more than four 
lensed images are detected. Based only on the position of the 
lensed images, an astrometric anomaly is detected as long as an 
astrometric accuracy smaller than 1 mas (but still compatible 
with the data error bars) is considered. Out of the last three 
systems (HE0435-1223, SDSS0924+0219, WFI2033-4723), 
constrained by relative astrometric positions derived from HST 
images, a good fit is obtained only for HE0435-12 23. A SIE+y 
mode l leads to x^r ~ 33 for SDSS0924-0219 (.Keeton et alJ 
I2OO6I) and to x^r ~ 15 for WFI2033-4723. In both cases, 
relaxing the constraint on the lens galaxy centroid allows to 
get a perfect fit of the astrometry. For WFI2033-4723 more 
complex models including a nearby group allow to reproduce 
the astrometry but they also need the light and ma ss distribution 
to be misahgned, which is not tot ally satisfactory (IVuissoz et alJ 
120081) . iKochanek & Dalall ( 12004 ^ compile seven quads (some 
of them also compiled here) for which they fitted SIE+y 
models. Unfortunately, there is only sparse information on the 
astrometric error used and we cannot infer any trend from this 
study. 

To conclude, out of the nine usable systems amongst 
the eleven quads, at least four show astrometric pertur- 
bations with respect to predictions of simple lens mod- 
els (B0128+437, MG0414+0534, SDSS0924+0219, WFI2033- 
4723). For three of the remaining systems (B 1608+656, 
MG2016+112, B2045+265), conclusions are diflicult to draw 
because the need to include a companion galaxy comes natu- 
rally from deep near-IR imaging. The last two systems (HE0435- 
1223, B1422+231) are well reproduced by simple smooth mod- 
els. 

Although the considered sample of quads gathers heteroge- 
neous data sets and analyses, it indicates that relative astrometry 
of quads often deviates from simple models expectation when 
trying to reproduce it to the mas precision. The considered sam- 
ple suggests that the situation is less critical for "central quads" 
(i.e. with the source lying close to the center of the central as- 
troid caustic) than for fold systems (i.e. source lying close to 
a fold caustic). This might be a normal geometrical effect (im- 
age positions vary more slowly when moving the source in the 
central region of the astroid caustic) but it might also be due 
to substructures lea ding to severe deformations of the caustics 
(iBradac et al.l 12002 ). It remains to be seen how significant this 
effect is with respect to the relative astrometric uncertainty on 
the image positions or the amount of shear. 



'" B 1933+503 also has VLA data but has already been ruled out, its 
configuration being easily reproduced with simple models because of 
the large uncertainty on the lens position. 

" The latter has been unveiled recently on AO images by 
lLagattutaet"^ ( l2010h . 



Substructures are not the only explanation for the frequent 
inability of simple lens models to fit the configuration of quads. 
Other possible explanations are astrometric perturbations due 
to the lens environment, asymmetries in the mass distribu- 
tion, disky/boxy projected mass profiles, offsets between the 
galaxy light centroid and mass c entroid. The last two so lutions 
seem however to be ruled out by Kochanek & Dalai (2004) and 
Yoo et al. (2006). The evidence that bright substructures/nearby 
satellite galaxies explain astrometric perturbations of some sys- 
tems suggest that substructures may be one of the major contrib- 
utors to the astrometric perturbations of quads. 

All this motivates a systematic study of the ability of simple 
models to reproduce the configuration of quads, with good con- 
trol on the error estimates and uniform modeling. Such a work is 
beyond the scope of this paper and is delayed to a forthcoming 
paper, when the iterative deconvolution method will have been 
applied to a larger number of quadruply imaged quasars. 



7. Conclusions 

In applying ISMCS, i.e. the MCS deconvolution algorithm com- 
bined with an iterative strategy, to HST/NIC2 images of seven 
lensed quasars, we have obtained accurate relative positional 
constraints on the lensed images, lensing galaxy as well as shape 
parameters for this galaxy. We reach an accuracy of around 1-2 
mas on the lensed image positions. We also detect for the first 
time a partial Einstein ring in two cases, the double HE 0047- 
1756 and the quadruple SDSS Jl 138+0314, and we highhght as 
well the already known ring in RX Jl 131-1231. Deeper images 
are needed to perform clear source reconstruction. In the case 
of HS 2209+1914, the deconvolved frame reveals a structure 
around the bulge of the lens galaxy which cannot be clearly iden- 
tified. This structure could be either an Einstein ring or some not 
well-resolved spiral arms of a big late-type galaxy. This question 
probably deserves further study: a spectrum of the suiTounding 
structures could give extra information about their true nature. 

We also obtain simple mass models for every sys- 
tem. In the case of doubles, both the isothermal and de 
Vaucouleurs profiles can reproduce the observed configuration. 
For SDSS Jl 155+6346, a good fit can only be reached with an 
extremely and anomalously high external shear, 0.392 in the case 
of a SIS and 0.453 in the case of a DV, indicating the presence 
of a galaxy group or cluster probably located outside the field of 
view of the ACS. For SDSS J1226-006, the large shear (r=0.1) 
is probably partially due to a nearby galaxy located 3'.'4 from the 
main deflector In the case of quads, a good;^'^ can only be ob- 
tained for one object: SDSS Jl 138+0314. The two other quads 
of our sample, RX J 1 1 3 1 - 1 23 1 and WFI J2026-4536, need more 
complicated models to account for their observed configuration. 
For RX Jl 131-1231, the offset between the light and mass dis- 
tribution cannot account for the astrometric perturbation we ob- 
serve. 

The study of the literature allows us to conclude that most 
of the quads cannot be modeled with simple profiles when the 
astrometic accuracy reaches around 1 mas: some need the pres- 
ence of companion galaxies, some others need substructures. 
This finding motivates the acquisition of mas astrometry for all 
the quads together with simple modeling. 

In the framework of the COSMOGRAIL collaboration, the 
next step for these seven systems is the acquisition of well- 
sampled light curves to extract time delays. Then, if the redshift 
of the lens is known, our astrometric constraints will help in re- 
ducing the systematic errors on the Hubble constant. 
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(a) HE 0047-1756 



■ 




(b) RX J 113 1-1231 



Fig. 1. HST/NIC2 original and deconvolved frames (resp. top left and top right), mean residual map from the first and from the last 
iteration of ISMCS (resp. bottom left and bottom right). 
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(d) SDSS J1I55+6346 



Fig. 1. continued. 
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(f) WFI J2026-4536 



m'^t-Ll.if.'rt^fJ^'J^-X 



Fig. 1. continued. 
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(g) HS 2209+1914 



Fig. 1. continued. 



